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ABSTRACT 

Successful  methods  for  numerically  processing  overpressure-time 
measurements  of  shock  and  blast  waves  and  sonic  booms  are  presented  and 
illustrated  with  sample  results.  These  methods  include:  a)  the  elimination 
of  high-frequency  random  and  other  noise  from  the  measurement  by  employing 
a  fast-Fourier-transform  technique  and  judiciously  altering  the  high- 
frequency  part  of  the  Fourier  spectrum;  b)  the  removal  of  low-frequency 
periodic  fluctuations  including  baseline  shift  and  drift  from  the  record  by 
using  specially  developed  techniques  described  in  the  report;  c)  the  extrac¬ 
tion  of  information  from  the  processed  shock-wave  record  including  shock- 
wave  arrival  time,  positive-phase  duration  and  impulse,  peak  overpressure, 
rise  time,  and  overpressure  decay  rate. 
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1.0  INTRODUCTION 

Good  quality  field  measurements  of,  for  example,  the  air  blast  wave 
from  a  gun  muzzle  or  a  chemical  or  nuclear  explosion,  the  shock  wave  from  a 
bursting  balloon  or  high-pressure  vessel,  and  the  sonic  boom  from  a  super¬ 
sonic  projectile  or  aircraft  are  normally  difficult  to  obtain,  because 
equipment  recording  conditions  in  the  field  are  generally  nonideal.  A  field 
overpressure  record  may  have  a  baseline  shift  and  also  drift.  Furthermore, 
the  record  may  contain  random  and  other  noise  from  instrumentation,  low- 
frequency  periodic  noise  from  power  lines,  extraneous  spikes  from  known  or 
unknown  origins,  and  additional  noise  from  other  sources  such  as  wind  turbu¬ 
lence  and  mechanical  vibration  of  the  gage  mount.  Even  when  the  most  modern 
instrumentation  is  used  in  a  field  trial,  the  past  problem  of  some  records 
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being  of  poor  quality  (i.e.,  noisy)  -is  still  with  us,  although  this  problem 
is  less  severe  today.  Past  measurements  of  shock  and  blast  waves  that 
exhibit  various  types  of  noise  including  baseline  shift  and  drift  can  be 
found,  for  example,  in  Refs.  1  to  6. 

A  poor  quality  field  measurement  cannot  always  bo  discarded  in 
favour  of  obtaining  a  new  and  better  one,  because  a  field  trial  cannot  al¬ 
ways  be  repeated  easily  or  economically.  Consequently,  each  less-than-ideal 
measurement  has  to  be  processed  carefully  to  eliminate  the  interference 
noise  or  its  effects  before  accurate  information  can  be  extracted  from  the 
record.  In  the  case  of  an  overpressure-versus-time  record  of  a  shock  or 
blast  wave,  this  information  could  include  the  shock-wave  arrival  time, 
positive-phase  duration  and  impulse,  peak  overpressure,  rise  time,  and  over¬ 
pressure  decay  rate. 

Ethridge's  method  (Ref.  2)  of  reading  and  smoothing  an  overpressure¬ 
time  record  of  a  blast  wave  is  a  combination  of  first  plotting  the  measured 
signature  on  semilogari thmic  graph  paper  (£n{Ap)  versus  t  and  Ap  versus  £n{tl) 
and  then  fitting  specially  selected  exponential  and  exponential -1 i ke  (Fried- 
lander)  functions  to  different  but  overlapping  segments  of  the  profile. 
Similarly,  Harvey  and  Fisher's  method  (Ref.  3)  of  smoothing  a  blast-wave 
measurement  consists  of  fitting  a  polynomial  equation  or  other  specially 
selected  function  (modified  Friedlander  expression)  to  the  entire  signature. 
For  the  case  of  a  blast-wave  record  containing  superposed  random  noise  and 
extraneous  spikes,  both  of  these  approaches  can  yield  fairly  good  estimates 
of  the  peak  overpressure,  decay  rate,  and  positive-phase  duration.  Further¬ 
more,  these  methods  give  a  smooth  analytical  expression  that  can  be  used 
easily  for  other  purposes,  such  as  in  a  calculation  of  the  positive-phase 
impulse. 

The  success  of  Ethridge1 '  and  Harvey  and  Fisher's  methods  depends 
to  a  large  extent  on  choosing  an  appropriate  function  to  represent  the 
measured  overpressure  signature.  In  the  case  of  a  blast  wave  from  a  conven¬ 
tional  chemical  or  nuclear  explosion,  an  exponential-like  function  is  known 
from  theory  and  experiment  to  be  a  good  choice.  For  the  more  general  case 
of  shock  and  blast  waves  and  sonic  booms  having  nonexponential  wave  forms, 
the  choice  of  a  suitable  function  to  yield  good  results  is  not  necessarily 


UNCLASSIFIED 


UNCLASSIFIED 


3 


straightforward.  Furthermore,  such  curve-fitting  methods  are  normally  not 
very  successful  when  the  overpressure  record  contains  baseline  drift  or  low- 
frequency  periodic  noise. 

A  method  of  processing,  evaluating,  and  extracting  air-blast 
information  from  an  overpressure  record,  which  does  not  utilize  curve¬ 
fitting  techniques,  has  been  developed  by  Vortman  and  Long  (Ref.  5).  In 
their  work,  baseline  shift  is  successfully  eliminated  by  subtracting  it  from 
the  record,  after  the  shift  or  offset  value  has  been  determined  from  a  con¬ 
sideration  of  both  the  overpressure  record  and  the  time  integral  of  this 
record.  Similarly,  linear  baseline  drift  is  eliminated  by  subtraction,  after 
the  rate  of  change  of  the  drift  has  been  determined.  High-frequency  noise  in 
a  record  having  an  unfavourable  signal-to-noise  ratio  is  removed  by  using  a 
specially  developed  low-pass  digital  filter.  A  special  peak-averaging  tech¬ 
nique  is  employed  which  effectively  removes  high-frequency  fluctuations  with¬ 
out  significantly  increasing  the  rise  time  of  the  steep-fronted  shock.  Low- 
frequency  periodic  noise  of  constant  amplitude  is  eliminated  from  the  record 
by  subtracting  a  suitable  periodic  signal  having  a  predetermined  amplitude 
and  frequency.  Blast-wave  information  was  then  retrieved  from  the  modified 
or  corrected  overpressure  signature.  Note  that  the  record  processing  and 
information  retrieval  are  done  automatically  with  a  digital  computer. 

Wells  (Ref.  7)  recently  investigated  the  usefulness  of  the  Fourier 
transform  in  processing  digital  data  to  remove  unwanted  frequency  content  or 
noise  from  a  measurement.  In  his  work  the  fast  Fourier  transform  is  used 
first  to  obtain  the  Fourier  amplitude  versus  frequency  spectrum,  from  which 
unwanted  frequency  content  can  be  removed  easily.  The  inverse  fast  Fourier 
transform  is  then  employed  to  obtain  the  modified  record.  It  is  shown  that 
a  signal  can  be  readily  isolated  from  both  high-frequency  noise  and  low- 
frequency  fluctuations  including  baseline  shift  and  drift,  provided  that  no 
significant  signal  frequencies  are  eliminated  with  the  noise.  In  other  words, 
the  success  of  this  method  depends  to  a  large  extent  on  how  well  the  dominant 
signal  frequencies  do  not  overlap  the  noise  frequencies. 

The  present  work  is  concerned  primarily  with  the  application  of 
signal-processing  methods  to  remove  baseline  shift  and  drift,  low-frequency 
periodic  noise,  and  high-frequency  random  and  other  noise  from  an  overpressure- 
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time  record  of  shock  and  blast  waves  and  sonic  booms.  Ideas  from  previous 
work  are  incorporated  freely  and  extended  when  necessary.  In  Chapter  2,  a 
method  which  successfully  removes  high-frequency  noise  from  a  measurement 
is  presented  and  illustrated  with  results.  This  method,  which  makes  use  of 
the  fast  Fourier  transform  and  its  inverse,  goes  beyond  the  past  technique 
of  removing  only  a  band  of  or  all  high-frequency  content  from  the  Fourier 
spectrum.  In  the  present  work  the  Fourier  amplitude  or  frequency  content  is 
judiciously  altered  so  that  the  high-frequency  content  associated  with  only 
the  noise  is  removed,  and  high-f requency  content  associated  with  the  over¬ 
pressure  signature  is,  therefore,  retained.  Then,  in  Chapter  3,  two  methods 
for  removing  baseline  shift  and  drift  and  low-frequency  periodic  noise  from 
a  record  are  presented  and  illustrated  with  results.  The  first  and  more 
general  method  involves  a  specialized  procedure  of  determining  the  baseline 
shift  and  the  amplitude,  frequency,  and  shape  of  the  periodic  noise  from  the 
baseline  signal  preceding  the  shock-wave  signature.  The  second  and  more 
limited  method  involves  fitting  a  special  function  consisting  of  a  constant 
(baseline  shift),  a  linear  time  term  (baseline  drift),  and  a  harmonic  term 
of  unknown  amplitude  and  phase  but  specified  frequency  (periodic  noise)  to 
the  baseline  signal  preceding  the  shock-wave  signature.  Then,  in  both  me¬ 
thods,  the  low-frequency  noise  can  be  subtracted  out  of  the  overpressure 
record. 

The  present  work  is  concerned  also  with  algorithms  for  retrieving 
information  from  the  modified  overpressure  measurement.  In  Chapter  4, 
methods  are  described  briefly  for  obtaining  the  shock-wave  arrival  time, 
positive-phase  duration  and  impulse,  peak  overpressure,  rise  time,  and  over¬ 
pressure  decay  rate. 

The  signal  processing  work  reported  herein  stems  mainly  from  having 
to  overcome  problems  encountered  in  the  computer  processing  of  numerous  sonic- 
boom  measurements  of  supersonic  projectiles  fired  from  military  guns 
(106-inp  M344A1 ,  3-inch/50  MK33,  3-inch/70  MK  34,  5-inch/S4  MK41).  These 
experimental  results  will  be  published  separately  as  DRES  reports.  Because 
of  nonideal  recording  conditions  encountered  during  the  field  trials  in  which 
overpressure  signatures  were  measured  at  different  distances  from  me  super¬ 
sonic  projectile,  some  of  the  overpressure  records  contained  various  types 
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of  noise.  Almost  all  of  the  records  had  some  degree  of  baseline  shift,  but 
only  a  few  had  noticeable  baseline  drift.  Although  very  few  of  the  records 
contained  significant  high-frequency  random  electronic  noise  with  an  un¬ 
favourable  signal -to-noise  ratio,  a  few  of  the  records  contained  significant 
high-frequency  mechanical  noise  from  gage-mount  vibration  or  ringing.  A  few 
other  records  had  marked  low-frequency  periodic  noise  of  about  60  Hz,  owing 
to  electrical  interference  from  110-volt  power  lines.  The  methods  presented 
in  Chapters  2  to  A,  to  remove  such  undesirable  noise  and  extract  accurate 
information  from  the  processed  records,  are  described  and  illustrated  with 
results  that  pertain  directly  to  sonic-boom  signatures  from  supersonic  pro¬ 
jectiles.  It  should  be  remembered,  however,  that  these  methods  are  much 
more  general  and  can  be  used  to  remove  similar  noise  and  extract  similar 
information  from  most  shock-  and  blast-wave  measurements. 

2.0  HIGH-FREQUENCY  filtering 

The  present  procedure  of  removing  high-frequency  noise  from  an 
overpressure-time  measurement  or  other  type  of  analog  record  consists,  essen¬ 
tially,  of  the  following  four  sequential  steps:  the  conversion  of  the  measure¬ 
ment  or  analog  record  into  a  digital  record;  the  use  of  the  Fourier  transform 
to  obtain  the  Fourier  amplitude  versus  frequency  spectrum;  the  judicious 
alteration  of  the  Fourier  spectrum  to  remove  undesirable  noise;  the  use  of 
the  inverse  Fourier  transform  to  obtain  the  modified  or  corrected  digital 
record  (free  of  undesirable  noise).  These  four  basic  steps  are  now  described 
in  more  detail . 

The  equipment  and  computer  programming  techniques  for  the  analog- 
to-digital  conversion  process  are  not  described  in  this  report,  because  they 
are  specific  to  the  type  of  tape  recorder  and  computer  available  to  the  user. 
It  is  worth  mentioning,  however,  that  the  sampling  rate  is  normally  chosen 
to  be  sufficiently  high  to  ensure  that  significant  high-frequency  content 
of  interest  is  retained  in  the  digital  record.  Other  considerations  such 
as  the  maximum  computer  storage  capability  or  computer  execution  cost  may 
sometimes  limit  the  maximum  sampling  rate. 

The  discrete  Fourier  transform  can  be  obtained  efficiently  from 
the  digital  record  by  using  a  relatively  new  algorithm  known  as  the  fast 
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Fourier  transform  (Refs.  8  to  10),  which  stems  from  the  original  work  of 
Cooley  and  Tukey  (Ref.  11).  A  brief  description  of  this  algorithm  and  the 
computer  program  subroutine  of  Ref.  12,  including  the  FORTRAN  subroutine 
listing,  are  given  for  interest  and  reference  in  Appendix  A. 

Undesirable  noise  in  the  discrete  Fourier  amplitude  versus 
freguency  spectrum  can  now  be  identified  and  removed.  For  the  simplest 
case  when  the  noise  contribution  to  the  Fourier  spectrum  is  distinct  from 
the  signal  contribution,  the  identification  and  removal  of  this  noise  con¬ 
tent  are  achieved  readily.  For  example,  consider  an  idealized  situation 
for  which  the  signal  contribution  to  the  Fourier  speccrum  occurs  at  low 
frequencies  and  the  noise  contribution  occurs  distinctly  at  higher  fre¬ 
quencies.  The  high-frequency  noise  content  can  then  be  removed  easily 
from  the  Fourier  spectrum  by  simply  setting  to  zero  all  of  the  noise  ampli¬ 
tudes  at  the  higher  frequencies.  In  the  converse  situation,  low-frequency 
noise  content  can  be  deleted  similarly  from  the  Fourier  spectrum  when  the 
signal  contribution  occurs  predominantly  at  higher  freauencies.  Furthermore, 
noise  content  occupying  a  distinct  frequency  band  or  bands  can  also  be 
deleted  easily. 

For  the  more  complex  case  when  the  noise  and  signal  contributions 
to  the  Fourier  spectrum  overlap  significantly,  identification  and  removal  of 
the  noise  may  then  be  quite  difficult  or  even  impossible.  If  all  of  the 
Fourier  amplitudes  in  a  frequency  band  that  contains  noise  content  were 
simply  set  equal  to  zero,  then  not  only  the  noise  contribution  but  also 
the  signal  contribution  in  this  band  would  be  deleted  from  the  Fourier 
spectrum.  In  order  to  circumvent  this  difficulty,  the  Fourie1'  amplitudes 
should  be  judiciously  altered  to  remove  only  the  noise  content  and,  there¬ 
fore,  retain  most  of  the  signal  content.  Such  a  Fourier  spectrum  alteration 
can  normally  be  achieved  only  if  the  noise  contribution  can  be  identified 
accurately.  Hence,  the  basic  shape  of  the  Fourier  spectrum  of  a  signal 
without  significant  noise  must  generally  be  known  beforehand. 

Noise  quite  often  causes  a  distinct  peak  or  peaks  in  an  otherwise 
fairly  smooth  Fourier  spectrum  of  a  signal.  If  the  basic  spectrum  shape  of 
the  noise-free  signal  is  known  in  advance,  these  noise  peaks  and  their  band- 
widths  are  then  recognized  easily.  Each  noise  peak  can  then  be  eliminated 
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from  the  frequency  domain  and  the  resulting  gap  bridged  by  a  smooth  poly¬ 
nomial  curve  which  is  joined  to  the  smooth  part  of  the  spectrum  on  each 
side  of  the  peak.  In  this  manner  the  dominant  noise  contribution  is 
effectively  removed,  while  the  signal  contribution  is  left  intact. 

The  last  of  the  four  steps  in  the  procedure  is  to  take  the  inverse 
Fourier  transform  of  the  altered  Fourier  spectrum  to  obtain  the  modified  or 
corrected  overpressure  record.  This  can  be  done  easily  by  using  the  same 
efficient  algorithm  for  the  fast  Fourier  transform,,  as  is  described  in 
Appendix  A. 

The  following  results  given  in  this  chapter  illustrate  the  use¬ 
fulness  of  the  foregoing  procedure  for  removinq  high-frequency  noise  from 
an  overpressure  measurement.  Firstly,  consider  an  ideal  sonic-boom  signa¬ 
ture  having  two  equal  amplitude  front  and  rear  shocks  with  instantaneous 
rise  times  and  separated  by  a  linear  decay  in  overpressure.  The  Fourier 
spectrum  as  calculated  by  the  fast  Fourier  transform  is  shown  in  Fig.  1, 
with  the  N-shaped  overpressure  signature  appearing  as  an  insert.  For  this 
and  all  other  results  presented  in  this  chapter,  the  analog-to-digital 
conversion  rate  is  96,000  points  per  second  and  the  Fourier  amplitudes  are 
calculated  at  512  equally  spaced  frequencies  between  zero  and  about  24  kHz. 

The  results  shown  in  Fig.  1  for  an  ideal  N-wave  of  similar  ampli¬ 
tude  and  duration  to  N-waves  measured  in  a  field  trial  *or  a  supersonic 
projectile  illustrate  some  of  the  basic  features  of  the  Fourier  spectrum. 

At  certain  frequencies  the  spectrum  is  peaked,  whereas  at  others  the  ampli¬ 
tude  drops  to  nearly  zero.  The  frequency  at  which  the  first  peak  occurs  is 
approximately  T"1  Hz  or  500  Hz,  as  might  be  expected  (Ref.  13),  where  T  is 
the  N-wave  duration.  Most  of  the  energy  spectral  density  of  the  N-wave  is 
contained  in  the  lower  frequencies,  less  than  about  10  T-1  Hz  or  5000  Hz. 

It  should  be  noted  that  analytical  results  are  available  in  Refs.  13,  14, 
and  15  for  the  Fourier  spectrum  or  energy  spectral  density  of  such  a  simple 
shaped  N-wave. 

The  Fourier  spectrum  for  a  measured  N-wave  from  a  supersonic 
projectile  (3-inch/50  MK33)  is  presented  in  Fig.  2.  The  original  over¬ 
pressure  signature  (insert)  differs  only  slightly  from  t  ,e  ideal  N-shape 
and  has  negligible  superposed  high-frequency  noise.  Consequently,  the 
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Fourier  spectrum  for  this  measured  signature  differs  only  to  a  small  extent 
from  that  for  an  ideal  N-wave  of  similar  amplitude  and  duration. 

In  Figure  3  the  Fourier  spectrum  of  a  measured  N-shaped  overpressure 
signature  (insert)  with  significant  high-frequency  noise  resulting  from  gage- 
mount  vibration  is  displayed.  By  comparing  this  spectrum  to  the  previous 
two  cases  for  the  ideal  and  measured  N-waves  (Figs.  1  and  2),  the  noise 
contribution  to  the  spectrum  can  be  identified  fairly  easily.  Additional 
energy  now  appears  in  the  frequency  bands  of  9  to  11  kHz  and  17  to  19  kHz. 
These  two  bands  correspond  directly  to  the  first  two  vibration  modes  of  the 
gage  mount. 

How  that  the  distortion  in  the  Fourier  spectrum  from  the  high- 
frequency  noise  has  been  identified,  it  can  be  removed  in  a  number  of  dif¬ 
ferent  ways.  The  simplest  metnod  would  be  to  set  the  Fourier  amplitude  to 
zero  throughout  each  of  the  affected  frequency  bands.  However,  in  this  case 
some  energy  not  associated  with  the  noise  would  also  be  deleted.  A  more 
sophisticated  approach  would  consist  of  using  a  polynomial  equation  to 
bridge  each  frequency  band,  being  matched  by  a  linear  least  squares  method 
to  the  spectrum  on  each  side  of  the  band.  In  the  present  work  it  was  found 
that  a  straight  line  was  quite  adequate  for  this  purpose. 

The  Fourier  spectrum  appearing  in  Fig.  3  has  been  reproduced  in 
Fig.  4,  with  the  gage-mount  ringing  removed  from  only  the  9  to  11  kHz  band. 
The  noise  in  the  17  to  19  kHz  band  is  not  very  significant.  After  applying 
the  inverse  fast  Fourier  transform  to  this  altered  spectrum,  one  obtains 
the  modified  overpressure  signature  shown  as  an  insert.  Much  of  the  high- 
frequency  ringing  no  longer  appears  in  this  modified  signature.  This  is 
more  clearly  illustrated  in  Fig.  5,  where  a  direct  comparison  is  made 
between  the  original  measurement  containing  gage-mount  ringing  and  the 
modified  record  with  the  first  mode  of  gage-mount  ringing  removed. 

Results  such  as  those  given  in  Fig.  5  aptly  demonstrate  that 
certain  high-frequency  noise  can  be  successfully  reduced  in  overpressure 
time  measurements  by  the  present  method.  The  success  of  this  method  de¬ 
pends  to  a  large  extent,  however,  on  how  accurately  the  noise  can  be 
identified  and  removed  from  the  Fourier  spectrum.  This  can  be  more  art 
than  science.  Noise  that  is  concentrated  in  a  narrow  frequency  band  or 
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bands  can  normally  be  deleted  most  easily. 

3.0  LOW-FREQUENCY  FILTERING 

Overpressure  measurements  from  field  trials  frequently  contain 
low-frequency  noise  including  baseline  shift  and  drift.  Baseline  shift 
and  drift  generally  stem  from  a  fixed  or  changing  voltage  bias  in  the 
recording  equipment  and  gage  response  to  changes  in,  for  example,  ambient 
temperature.  Periodic  noise  can  be  the  result  of  interference  from  110- 
volt  power  lines  or  other  electronic  wiring  and  equipment.  In  the  present 
work  the  periodic  noise  had  a  fixed  frequency  of  about  60  Hz  and  was  sinu¬ 
soidal  in  shape.  However,  periodic  noise  having  a  nonsinusoidal  shape  is 
also  very  often  encountered  in  practice. 

Two  methods  were  devised  to  remove  low-frequency  noise  from  an 
overpressure  measurement.  The  first  method  is  quite  general  and  has  the 
capability  of  numerically  filtering  out  of  a  digitized  record  both  base¬ 
line  shift  and  low-frequency  periodic  noise  having  almost  any  unknown  shape 
and  frequency.  The  second  method  is  simpler,  requires  significantly  less 
computer  storage  and  execution  time,  but  is  more  limited  in  capability. 

Both  baseline  shift  and  drift  and  also  low-frequency  sinusoidal  shaped 
noise  of  specified  frequency  (e.g.,  60  Hz)  can  be  numerically  filtered 
out  of  an  overpressure  record.  These  two  methods,  labelled  A  and  B,  are 
described  and  illustrated  with  graphical  results  in  the  following  two 
sections. 

3.1  Method  A 

Consider  an  overpressure  record  from  a  field  trial,  which  has  an 
initial  length  of  baseline  signal  preceding  the  shock-wave  signature 
(e.g.,  see  Fig.  8a).  Suppose  that  this  baseline  signal  ahead  of  the  shock- 
wave  signature  contains  both  a  baseline  shift  and  low-frequency  noise  which 
is  nonsinusoidal  but  periodic.  This  shift  and  low-frequency  noise  is 
assumed,  of  course,  to  continue  through  the  shock-wave  signature  to  the  end 
of  the  record.  Method  A  consists,  essentially,  of  using  the  first  part  of 
the  overpressure  measurement  preceding  the  shock-wave  signature  to  estab¬ 
lish  one  cycle  of  a  reference  baseline  signal  that  accurately  represents 
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the  baseline  shift  and  low-frequency  periodic  noise.  The  resulting  refe¬ 
rence  signal  can  then  be  subtracted  from  the  overpressure  measurement,  in 
order  to  remove  the  baseline  shift  and  periodic  noise,  leaving  a  modified 
overpressure  record  containing  a  much  more  definitive  shock-wave  signature 
with  only  superposed  random  noise  (e.g.,  see  Fig.  8b). 

The  first  step  in  Method  A  is  to  determine  the  number  of  points 
in  the  discrete  overpressure  record  or  baseline  signal  preceding  the  shock- 
wave  signature.  This  number  which  is  denoted  by  the  symbol  Mp  can  be 
obtained  fairly  easily  by  conducting  a  computer  search  for  the  arrival  of 
the  shock  wave  in  the  overpressure  record. 

The  second  step  is  to  establish  an  initial  reference  signal  to 
represent  the  shift  and  periodic  noise  in  the  baseline  signal.  The  periodic 
noise  frequency  is  first  guessed  or  estimated,  making  certain  that  the  cho¬ 
sen  value  is  overestimated  slightly  (e.g.,  65  Hz  if  the  noise  frequency  is 
about  60  Hz).  For  an  overpressure  record  stored  in  discrete  form  (e.g., 
96,000  points  per  second),  the  total  number  of  points  lip  in  one  cycle  of 
this  initial  reference  signal  can  then  be  obtained  (e.g.,  1476  points  for 
a  frequency  guess  of  65  Hz).  The  successive  points  from  the  start  of  the 
overpressure  record  (N  =  1)  to  the  No-th  point  (N  =  Nq),  therefore, 
establishes  the  first  cycle  of  the  initial  reference  signal.  Successive 
duplication  of  these  values  then  establishes  subsequent  cycles  of  the 
initial  reference  signal,  having  a  total  number  of  points  equal  to  Np  or 
greater. 

The  third  step  is  to  compare  the  first  Np  values  of  the  reference 
signal  to  the  first  Np  values  of  the  overpressure  record.  A  convenient 
indicator  of  how  well  these  two  signals  agree  is  the  well-known  standard 
deviation: 

b 

°o  =  l  n^/Nb.  (1) 

N=1 

If  each  N-  h  difference  between  the  reference  signal  and  overpressure 
record  was  zero,  then  the  standard  deviation  would  also  be  equal  to  zero, 
and  the  correlation  between  the  two  signals  would  be  perfect.  Since  the 
initial  frequency  guess  was  slightly  high,  however,  the  standard  deviation 
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in  this  first  case  would  be  nonzero. 

The  fourth  step  is  to  pick  a  slightly  lower  frequency  than  the 
first  estimate,  and  then  repeat  the  comparison  to  obtain  a  new  standard 
deviation.  This  is  achieved  easily  by  using  Nq  +  1  or  N-j  successive 
values  of  the  overpressure  record  as  one  cycle  of  a  new  reference  signal. 

For  this  case  with  N-j  points  per  cycle,  let  the  new  value  of  the  standard 
deviation  be  denoted  by  oj. 

The  fifth  step  is  to  continue  the  procedure  of  slowly  increasing 
the  number  of  points  in  one  cycle  of  the  reference  signal  (or  decreasing 
the  frequency)  and  monitoring  the  new  values  of  the  standard  deviation. 

As  the  frequency  of  the  reference  signal  decreases  and  approaches  the 
noise  frequency,  the  standard  deviation  diminishes  slowly  at  first,  then 
quite  rapialy,  and  finally  levels  out  at  small  nonzero  values  (see  Fig.  6). 
When  the  frequency  of  the  reference  signal  is  decreased  further,  below  the 
noise  frequency,  the  standard  deviation  begins  eventually  to  increase 
rapidly  from  the  small  nonzero  values,  and  then  increases  more  slowly.  The 
fairly  flat  minimum  in  the  resulting  standard-deviation  distribution  is 
nonzero  (Fig.  6),  because  perfect  matching  of  the  reference  signal  and  over¬ 
pressure  record  is  unachievable,  owing  to  random  noise  in  the  overpressure 
record.  Furthermore,  such  random  noise  also  causes  small  erratic  variations 
in  the  discrete  standard-deviation  values  in  the  vicinity  of  the  minimum. 
Although  the  exact  minimum  is  difficult  to  locate  precisely,  a  good  estima¬ 
tion  can  be  made.  For  example,  a  parabolic  curve  can  be  fitted  by  a  linear 
least-squares  method  to  the  standard-deviation  values  in  the  vicinity  of 
the  minimum  (Fig.  6),  and  the  location  of  the  minimum  can  then  be  obtained 
from  this  parabola.  This  location  yields  directly,  of  course,  a  good 
estimate  of  the  frequency  of  the  periodic  noise. 

The  sixth  step  is  to  refine  the  previous  frequency  estimate  of 
the  periodic  noise.  A  straight  line  is  first  fitted  by  a  linear  least- 
squares  method  to  a  few  values  at  the  beginning  of  the  reference  signal. 

This  straight  line  is  then  used  to  yield  an  average  amplitude  at  the  start 
of  the  reference  signal.  The  effects  of  random  noise  on  this  average  ampli¬ 
tude  are,  therefore,  minimized.  Similarly,  a  straight  line  is  fitted  to  a 
few  values  at  the  end  of  the  first  cycle  of  the  reference  signal,  and  an 
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average  amplitude  is  also  obtained  at  this  location.  If  these  two  ampli¬ 
tudes  at  the  start  and  end  of  the  first  cycle  are  found  to  be  different, 
the  following  adjustment  to  the  frequency  is  then  made.  The  end  of  the 
first  cycle  of  the  reference  signal  is  lengthened  or  shortened  such  that 
the  average  amplitude  at  the  end  of  the  first  cycle  changes  to  match  exact¬ 
ly  the  average  amplitude  at  the  beginning.  Each  cycle  of  the  reference 
signal,  therefore,  starts  and  ends  with  the  same  average  amplitude.  The 
preceding  adjustment  to  the  periodic  noise  frequency  is  normally  quite 
small,  but  it  can  be  an  important  step  in  determining  a  precise  value  of 
the  frequency. 

The  seventh  step  is  to  obtain  a  final  reference  signal  for  the 
periodic  noise,  having  the  previously  determined  frequency.  A  number  of 
sequential  cycles  of  periodic  noise  in  the  overpressure  record,  starting 
with  the  first  cycle,  are  now  averaged  to  obtain  a  final  reference  signal. 
The  number  of  cycles  that  is  employed  is  made  as  large  as  possible,  utili¬ 
zing  almost  all  of  the  available  baseline  signal  preceding  the  shock-wave 
signature.  Such  an  averaging  procedure  helps  to  minimize  the  effects  of 
anomalous  noise  (e.g.,  an  infrequent  spike)  that  might  have  occurred  in 
the  first  cycle  of  the  overpressure  record. 

The  eighth  and  final  step  is  to  now  use  the  final  reference  signal 
to  remove  the  periodic  noise  from  the  overpressure  record.  This  is  done 
quite  simply  by  subtracting  the  reference  signal  from  the  overpressure 
record,  leaving  a  modified  overpressure  record  containing  an  insignificant 
baseline  shift  and  a  relatively  unfluctuating  baseline  signal  with  a  well- 
defined  shock-wave  signature,  along  with  some  minor  superposed  random  noise. 

The  capability  of  Method  A  is  well  demonstrated  b>  the  two  sets  of 
results  given  in  Figs.  7  and  8.  In  each  set  the  original  overpressure 
measurement  appears  at  the  top  of  the  figure  and  the  modified  record  is 
shown  below.  From  these  two  sets  of  results  it  can  be  clearly  seen  that 
Method  A  is  very  successful  in  removing  the  baseline  shift  and  undesirable 
low-frequency  periodic  noise. 

For  the  sake  of  interest  the  overpressure  measurement  shown  in 
Fig.  7  is  a  shock-wave  record  from  a  supersonic  projectile  (3- inch/50  MK33). 
Its  speed  was  460  m/s,  corresponding  to  a  Mach  number  of  1.4,  and  the 
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measurement  was  taken  at  a  distance  of  15  m  to  one  side  of  the  flight  path. 
The  shock-wave  overpressure  measurement  shown  in  Fig.  8  shows  only  the 
first  part  of  the  muzzle  blast  wave  from  the  artillery  gun  that  fired  the 
supersonic  projectile.  This  measurement  was  taken  185  m  out  front  of  the 
gun  muzzle. 

3.2  Method  B 

Method  B  can  be  used  to  eliminate  baseline  shift,  drift  and  low- 
frequency  periodic  noise  from  an  overpressure  measurement.  The  periodic 
noise  is  assumed  to  be  sinusoidal  in  shape  with  a  known  frequency.  The 
following  equation  is  used  to  represent  the  low-frequency  noise  including 
baseline  shift  and  drift. 

N(t)  =  a  sin(2iTft  +  $)  +  f*  +  yt  (2) 

The  constant  a  is  an  initially  unknown  amplitude  of  the  periodic  noise, 
f  denotes  the  known  frequency  of  the  periodic  noise,  t  denotes  time,  and 
4>  is  an  initially  unknown  phase  angle.  The  constant  fi  denotes  an  initially 
unknown  value  of  the  baseline  shift,  and  the  constant  y  is  an  initially 
unknown  rate  of  change  of  baseline  drift.  Baseline  drift,  therefore,  is 
assumed  to  be  linear. 

Equation  2  can  be  expressed  in  the  following  alternate  form,  as  a 
combination  of  linear  functions,  which  is  required  for  the  present  work. 

N(t)  =  a  cos ( <t>)  sin(2nft)  +  u  s i n( <j>)  cos(2nft)  +  e  +  yt  (3) 

A  linear  least-squares  method  can  then  be  employed  to  fit  this  equation  to 
the  discrete  baseline  data  preceding  the  shock-wave  signature,  resulting 
in  a  direct  determination  of  the  following  four  coefficients:  a  cos(<j>), 
a  sin(<t>),  (3,  and  y.  From  the  values  of  the  first  two  coefficients,  indi¬ 
vidual  values  of  a  and  $  can  now  be  obtained. 

Now  that  the  constants  a,  6,  y,  and  $  have  been  determined,  Eq.  2 
or  3  can  be  used  to  eliminate  the  sinusoidal  noise  and  baseline  shift  and 
drift  from  the  overpressure  record.  This  is  accomplished  easily  by  sub¬ 
tracting  the  noise  signal  (Eq.  2  or  3)  from  the  overpressure  record.  The 
resulting  modified  record  should  have  insignificant  baseline  shift  and 
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drift  and  consist  of  a  relatively  unfluctuating  baseline  signal  with  a  well- 
defined  shock-wave  signature. 

The  capability  of  Method  B  is  aptly  demonstrated  by  the  two  sets  of 
results  given  in  Figs.  9  and  10.  The  original  set  of  overpressure  measure¬ 
ments  appears  first  in  Fig.  9.  In  these  measurements  the  periodic  noise  is 
most  severe  in  the  second  record,  whereas  the  worst  baseline  shift  occurs 
in  the  third  record.  The  set  of  modified  overpressure  records  is  shown  for 
comparison  in  Fig.  10.  It  can  be  readily  seen  that  the  periodic  noise 
(assumed  to  be  60  Hz)  and  baseline  shift  have  been  almost  entirely  removed. 
Some  minor  low-frequency  noise  remains  superposed  on  only  the  second  record. 
In  conclusion.  Method  B  is  quite  effective  in  removing  such  noise  from  an 
overpressure  measurement. 

In  the  process  of  fitting  Eq.  2  to  the  discrete  overpressure  data, 
in  order  to  determine  the  values  of  a,  6,  y,  and  $,  data  for  a  short  time 
interval  of  only  5  ms  was  used.  This  corresponds  to  only  three-tenths  of 
a  cycle  of  the  periodic  noise.  Such  a  small  tvne  interval  was  found  to  be 
quite  satisfactory  for  our  overpressure  measurements  for  supersonic  pro¬ 
jectiles.  A  larger  time  interval,  when  available,  produced  only  marginal 
improvements  in  removing  the  periodic  noise  and  baseline  shift.  Further¬ 
more,  it  was  discovered  that  the  baseline  drift  (yt)  relative  to  the  base¬ 
line  shift  (b)  was  negligible  over  the  short  time  interval  of  interest 
(10  to  40  ms).  As  a  consequence,  only  the  baseline  shift  (b)  needed  to  be 
taken  into  account,  and  the  value  of  the  rate  of  change  of  baseline  drift 
(y)  could  be  set  equal  to  zero  from  the  outset. 

Methods  A  and  B  should  be  considered  to  be  complementary.  When 
only  one  cycle  or  a  fraction  of  a  cycle  of  periodic  noise  preceding  the 
shock-wave  signature  is  available  for  defining  the  reference  signal,  then 
Method  B  has  to  be  employed.  When  two  or  preferably  more  cycles  of  noise 
are  available,  Metnod  A  should  be  used.  Since  Method  A  is  more  general  in 
that  nonsinusoidal  noise  with  an  unknown  frequency  can  be  handled,  it  is 
normally  superior  in  eliminating  low-frequency  periodic  noise. 

For  the  sake  of  interest,  the  overpressure  measurements  presented 
in  Fig.  9  are  of  the  shock  wave  from  a  supersonic  projectile  (3-inch/ 

70  MK34).  This  projectile  was  moving  3  m  above  and  parallel  to  the  ground 
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at  a  speed  of  610  m/s,  corresponding  to  a  Mach  number  of  1.8.  The  four 
measurements  labelled  a,  b,  c,  and  d  were  recorded  at  different  distances 
of  7,  15,  36,  and  42  m,  respectively,  to  one  side  of  the  flight  path.  In 
the  top  three  records  of  Figs.  9  and  10  the  N-shaped  shock  wave  is  followed 
by  another  N-wave.  This  latter  N-wave  is  a  direct  consequence  of  the 
ground  acting  as  a  reflecting  plane.  In  the  bottom  record,  the  initial 
and  reflected  N-waves  just  begin  to  overlap. 

4.0  INFORMATION  RETRIEVAL  FROM  AN  OVERPRESSURE  RECORD 

Numerous  overpressure  records  and  other  signals  very  often  need 
to  be  processed  quickly  and  efficiently  to  obtain  various  types  of  infor¬ 
mation.  Automated  computer  processing  and  extraction  of  information  is 
most  often  the  only  feasible  approach.  Although  this  approach  is  usually 
quick  and  convenient  for  a  computer-program  user,  precise  definitions  and 
algorithms  must  be  developed  initially  to  ensure  that  the  user  obtains 
accurate  numerical  information  from  the  records.  In  the  following  six 
sections  of  this  chapter  our  methods  of  obtaining  information  from  over¬ 
pressure  records  for  shock-wave  arrival  tine,  positive-phase  duration  and 
impulse,  peak  overpressure,  rise  time,  and  initial  overpressure  decay  rate 
are  described  briefly.  The  methods  are  not  necessarily  new;  however,  they 
are  presented  herein  for  both  completeness  and  future  reference.  Research¬ 
ers  in  the  past  have  not  documented  the  details  •'  their  methods. 

4.1  Time  of  Arrival 

The  time  of  arrival  of  the  shock  wave  relative  to  a  reference  time 
is  determined,  essentially,  in  the  following  manner.  Near  the  beginning 
of  each  overpressure  record  an  initial  timing  mark  is  normally  available 
as  a  convenient  reference  point  (e.g.,  see  Fig.  11).  Starting  from  this 
reference  point  (time  tg),  or  the  beginning  of  the  record  in  the  absence 
of  such  a  mark,  the  discrete  overpressure  data  can  be  examined  numerically, 
in  a  sequential  point-by-point  manner.  Eventually,  a  few  consecutive  over¬ 
pressure  values  will  consistently  increase  (at  time  tj) ,  thereby  marking 
the  first  few  points  of  the  front  shock  of  the  N-wave.  The  difference 
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between  the  times  corresponding  to  the  first  value  of  the  ascending  sequence 
and  the  reference  point  (tj-tg)  is,  therefore,  the  time  of  arrival  of  the 
shock  wave. 

Although  this  procedure  is  simple  and  correct  conceptually,  addi¬ 
tional  computational  steps  or  checks  are  required  in  practice  to  avoid 
erroneous  results.  Firstly,  low-level  random  noise  in  the  overpressure 
record  may  contain  a  sufficiently  long  run  of  ascending  overpressure  values 
(Fig.  11),  which  could  stop  the  point-by-point  search  prematurely  (at  t^ ) . 

In  order  to  circumvent  this  problem  a  threshold  overpressure  level  can  be 
specified.  The  search  then  stops  only  when  the  run  of  ascending  values 
all  exceed  this  threshold.  Secondly,  an  infrequent  spike  in  an  overpres¬ 
sure  record  can  also  cause  a  premature  stopping  (at  t^),  yielding  an 
erroneous  time  of  arrival.  To  avoid  this  problem  the  shape  of  the  following 
signature  can  be  explored.  The  point-by-point  search  is  stopped  only  if  the 
signature  is  N-shaped.  Other  additional  safeguards  to  prevent  erroneous 
results  sometimes  have  to  be  devised  to  overcome  problems  peculiar  to  a 
particular  record. 

4.2  Positive-phase  Duration 

The  duration  of  the  initial  positive  phase  of  the  overpressure 
signature  (Fig.  11)  is  defined  as  the  time  interval  from  the  arrival  of 
the  N-wave  (t^)  to  the  time  (tg)  when  the  overpressure  signature  crosses 
the  horizontal  axis  and  becomes  negative  (i.e.,  tg-t^l.  The  latter  time 
(tg)  must  be  chosen  carefully  from  the  discrete  overpressure  data  so  that 
extraneous  crossings  (e.g.,  at  t^)  resulting  from  random  or  other  noise 
are  excluded.  One  common  procedure  of  obtaining  the  crossing  time  accu¬ 
rately  is  to  first  fit  a  straight  line  to  the  data  in  the  vicinity  of  the 
crossing.  Then  the  crossing  time  can  be  obtained  from  this  straight  line 
with  confidence.  Another  common  procedure  is  to  first  obtain  the  time 
integral  of  the  overpressure  record,  starting  the  integration  from  the 
arrival  time  (tj).  The  crossing  time  then  corresponds  to  the  time  at  which 
the  integral  is  a  maximum.  These  two  procedures  normally  yield  very  simi¬ 
lar  results  for  the  positive-phase  duration.  In  the  authors'  work  with 
shock  waves  from  supersonic  projectiles  the  latter  integral  procedure  is 
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employed. 


4.3  Positive-phase  Impulse 


The  positive-phase  impulse  (I)  is  defined  as  the  value  of  the  tine 
integral  of  the  positive-phase  overpressure  (Ap)  from  the  N-wave  arrival 
time  (t^)  to  the  end  of  the  positive-phase  duration  (t^),  as  shown  below. 


I 


tg 

.  /  Ap(t)  dt 

z3 


(4) 


The  integration  of  the  discrete  overpressure  data  by  using  a  trapezoidal 
or  Simpson's  rule  is  quite  straightforward.  Note  that  the  impulse  value 
is  normally  obtained  simultaneously  with  the  determination  of  the  end  of 
the  positive-phase  duration  (tg)  when  the  integral  method  is  used. 

4.4  Peak  Overpressure 

The  peak  overpressure  of  the  front  shock  of  an  N-wave  is  easily 
picked  from  the  discrete  overpressure  data.  The  time  at  which  this  peak 
occurs  is  also  obtained  easily  (t^  in  Fig.  11).  To  minimize  the  effects 
of  random  noise  on  this  peak  value,  it  can  be  modified  slightly  by 
averaging  it  with  a  few  preceding  and  following  overpressures .  This 
averaging  procedure  is  meaningful  only  when  the  digitization  rate  is 
sufficiently  high  so  that  the  averaged  overpressures  are  all  an  essential 
part  of  the  peak.  The  peak  overpressure  obtained  in  this  mariner  is  denoted 
by  the  symbol  aPq  for  future  reference. 

An  alternate  estimate  of  the  peak  overpressure  can  be  obtained  in 
the  following  wa;, .  An  exponential  curve  is  first  fitted  to  the  discrete 
overpressure  data  in  the  time  interval  starting  at  t^  and  ending  at 
(t4  +  tg)/2  (Fig.  11).  The  peak  overpressure  is  then  obtained  directly 
from  this  exponential  curve  at  the  time  t^,  and  denoted  by  the  symbol  Ap^ . 
In  the  author's  work  with  shock  waves  from  supersonic  projectiles,  both  of 
the  peak  overpressures  (Ap^  and  Ap-j )  were  determined  and  found  to  be  in 
good  agreement. 
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4.5  Rise  Time 

The  rise  time  of  the  front  shock  of  the  N-wave  is  defined  as  the 
time  interval  starting  from  the  arrival  of  the  N-wave  (t^)  and  ending  at 
the  time  (t^)  when  the  overpressure  is  a  maximum  (i.e.,  t^-t^).  From  the 
previously  determined  data  for  t^  and  the  rise  time  can  be  obtained 
easily. 

Other  definitions  of  shock-wave-  rise  time  are  encountered  in  the 
open  literature.  For  the  fully  dispersed  shock  structure  of  weak  waves 
and  impulse  sound  effects  on  hearing  the  rise  time  is  often  defined  as 
1.25  times  the  time  for  the  overpressure  to  rise  from  10  to  90%  of  its 
peak  value  (Ref.  16).  The  rise  time  obtained  with  this  definition  is 
generally  10  to  30%  shorter  than  in  the  present  case. 

4.6  Overpressure  Decay  Rate 

The  decay  rate  of  the  overpressure  just  behind  the  front  shock  of 
the  N-wave  is  determined  directly  from  the  exponential  curve  which  was  fit 
ted  previously  to  the  discrete  overpressure  data  in  Section  4.4.  This 
decay  rate  is  equal  to  the  derivative  of  the  overpressure  at  the  time  t^ 
where  the  peak  overpressure  ap-|  was  taken. 

4.7  Sample  Results 

The  previous  methods  of  obtaining  the  shock-wave  time  of  arrival, 
positive-phase  duration  and  impulse,  peak  overpressure,  rise  time,  and 
overpressure  decay  rate  have  been  applied  to  the  four  shock-wave  over¬ 
pressure  signatures  shown  in  Fig.  10.  The  results  are  summarized  for 
interest  in  Table  1.  Note  that  the  information  has  been  extracted  from 
the  modified  overpressure  records  from  which  the  baseline  shift  and  low- 
frequency  and  high-frequency  noise  have  been  removed.  This  ensures  that 
the  information  is  of  the  highest  possible  degree  of  accuracy. 

5.0  CONCLUDING  REMARKS 

Numerical  methods  of  successfully  processing  digitized  overpres- 
suro-versus-time  measurements  of  shock  waves  from  supersonic  projectiles 
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have  been  presented,  illustrated  with  results,  and  discussed  in  detail. 

It  has  been  shown  that,  once  hinh-frequency  noise  is  identified  in  the 
Fourier  spectrum  of  the  measurement,  this  noise  can  then  be  removed 
successfully  by  judiciously  altering  the  spectrum  before  applying  the 
inverse  Fourier  transform.  As  a  prerequisite  to  the  spectrum  alteration, 
a  fair  knowledge  of  the  spectrum  shape  without  noise  effects  is  normally 
required. 

It  has  also  been  shown  that  baseline  shift  and  drift  and  low- 
frequency  periodic  noise  can  be  removed  successfully  from  a  measurement  by 
two  complementary  methods.  In  Methods  A  and  B  baseline  shift  is  eliminated 
equally  well,  but  only  Method  B  has  the  capability  of  removing  baseline 
drift.  Method  A  is  generally  superior  to  Method  B  for  removing  low- 
frequency  periodic  noise,  since  it  was  devised  to  remove  periodic  noise  of 
initially  unknown  shape  and  frequency.  Method  B  was  devised  to  remove 
periodic  noise  w;th  a  known  sinusoidal  shape  and  specified  frequency.  It 
should  be  remembered  that,  in  order  to  use  Method  A,  two  or  preferably  more 
cycles  of  periodic  noise  in  the  baseline  signal  preceding  the  shock-wave 
signature  are  required  for  characterizing  the  noise  shape  and  frequency. 

When  only  one  cycle  or  a  fraction  of  a  cycle  of  periodic  noise  is  available, 
then  Method  B  has  to  be  employed. 

The  present  methods  of  removing  low-  and  high-frequency  noise  and 
retrieving  accurate  information  from  a  modified  overpressure  measurement 
have  been  applied  in  the  data  reduction  of  numerous  overpressure  measure¬ 
ments  of  the  shock  wave  from  supersonic  projectiles  fixed  from  military  guns 
(106-mm  M344A1,  3-inch/50  MK33,  3-i nch/70  MK34,  5-inch/54  MK41).  These 
experimental  results  will  be  published  separately  as  DRES  reports.  It  is 
important  to  note  that  these  methods  are  not  restricted  in  use  to  only 
shock-wave  records  from  supersonic  projectiles,  but  can  also  be  applied  to 
other  similar  records.  For  example,  the  methods  are  almost  directly  appli¬ 
cable  to  shock-wave  measurements  from  both  shock  tubes  and  blast-wave 
simulators  and  also  to  blast-wave  records  from  chemical  explosions. 
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TABLE  I 


DATA  FROM  A  SHOCK  WAVE  OF  A  SUPERSONIC  PROJECTILE 
( 3- i nch/70  MK34,  projectile  speed 
and  Mach  number  are  610  m/s  and  1.0) 


Shock-Wave  Properties 

Shock-Wave  Data  for  Signatures  shown  in: 

Fig.  10a 

Fig.  10b 

[ 

Fig.  10c 

Fig.  lOd 

I  illative  time  of 
arrival  (ms) 

0 

20 

72 

88 

Positive-phase  duration 
(ms) 

0.698 

0.823 

1.04 

1.08 

Positive-phase  impulse 
(kPa  ms) 

0.554 

0.377 

0.222 

0.192 

Peak  overpressure 
ap0  ( kPa) 

1 .77 

0.960 

0.398 

0.344 

apt  (kPa) 

1.71 

0.955 

0.406 

0.336 

Shock  rise  time  (us) 

31 

41 

72 

72 

Initial  decay  rate 
(kPa/ms) 

-4.10 

-1.85 

-0.495 

-0.405 

UNCI  ASS  IF!  ED 
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FIG.  I.  FOURIER  SPECTRUM  OF  AN  IDEAL  N-WAVE. 
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AMPLITUDE  (Pcms) 


FREQUENCY  (kHz) 


FIG.  2.  FOURIER  SPECTRUM  OF  AN  N-WAVE 
FROM  A  SUPERSONIC  PROJECTILE. 
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FREQUENCY  (kHz) 


FIG.  3.  FOURIER  SPECTRUM  OF  AN  N-WAVE  WITH 
SUPERPOSED  GAGE -MOUNT  RINGING. 
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FIG.  4.  FOURIER  SPECTRUM  OF  AN  N-WAVE  WITH 
THE  FIRST  MODE  OF  THE  SUPERPOSED 
GAGE-MOUNT  RINGING  REMOVED. 
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FIG.  9 


ORIGINAL  OVERPRESSURE  RECORDS  OF  THE 
SHOCK  WAVE  FROM  A  SUPERSONIC  PROJECTILE 
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MODIFIED  OVERPRESSURE  RECORDS  OF  THE 
SHOCK  WAVE  FROM  A  SUPERSONIC  PROJECTILE. 
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FIG.  II.  SAMPLE  SHOCK-WAVE  OVERPRESSURE  RECORD. 
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APPENDIX  A 

FAST  FOURIER  TRANSFORM  SUBROUTINE 


A  fast-Fourier-transform  subroutine,  similar  to  that  in  Ref.  A1 ,  is 
used  to  compute  the  discrete  Fourier  transform  by  the  nlogn  method.  This 
subroutine  is  called  FFT,  and  a  sample  FORTRAN  listing  is  given  at  the  end 
of  this  appendix.  The  extremely  efficient  nlogn  method  requires  that  the 
data  consists  of  q  values,  where  q  equals  2n  and  n  is  some  positive  integer. 
Subroutine  FFT  does  not  make  use  of  any  other  subroutines.  Its  call  state¬ 
ment  is  simply 

CALL  FFT(N,  X,  SIGN). 

The  subroutine  inputs  are 

N  =  n, 

X  =  Xq,  Xr  X2,  . .  Xq_.|  (complex  numbers) 

SIGN  --  -1.0  or  +1.0. 


The  subroutine  output  is 

X  =  Yq,  Y.j,  Y^,  . . Y^  (complex  numbers), 

where  , 

q- 1 

Y.  -  l  X.  exp(2nikj/q) ,  if  SIGN  =  -1.0, 

J  k-0 


Y.  =  1  l  X.  exp(-2uikj/q),  if  SIGN  =  +1.0. 

J  q  k=0 

If  X^  (k  «  0,  1,  ...  q-1)  represents  a  time  function,  then  the  sub¬ 
routine  output  Y^  (k  =  0,  1,  ...  q-1)  in  the  case  of  SIGN  -  -1.0  represents 
the  frequency  spectrum  at  frequencies  f  =  j/q  cycles  per  time  unit,  where 
j  =  0,  1 ,  . . .  q-1  .  In  particular,  if  X^  is  real  valued,  then 

q-1 

Re  { Y . }  =  l  X.  cos(2irkj/q) 

J  k=0 


is  the  cosine  spectrum,  and 
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q-i 

Im{Y .}  =  -  l  X.  sin(2itkj/q) 

J  k=0 

is  the  sine  spectrum.  To  recover  the  time  function  X^  we  may  call  Subroutine 
FFT  again,  but  this  time  with  SIGN  =  +1.0,  That  is,  given  the  array  X,  the 
call  statements 

CALL  FFT(N,  X,  -1.0), 

CALL  FFT ( N ,  X,  +1.0), 

one  after  another,  yield  the  same  X  as  final  output,  because 
q-l  [  q-l 

i  I  l  X.  exp(-2irikj/q)  exp(2ni jm/q)  =  X  . 
q  j-o  [  k=0  K  m 

Likewise,  given  the  array  X,  the  call  statements 
CALL  FFT(N,  X,  +1.0), 

CALL  FFT(N,  X,  -1.0), 

one  after  another,  yield  the  same  X  as  final  output. 

As  a  numerical  example,  suppose  we  have  the  function 

X  =  (1,  2,  3,  A). 

Then  the  call  statement 

CALL  FFT(2,  X,  +1.0) 

would  yield  the  Fourier  transform 

X  =  (2.5,  -0.5-0. 5.,  -0.5,  0.5+0.51). 

If  we  execute  the  call  statement 

CALL  FFT(2,  X,  -1.0), 

we  return  to  the  given  function 

X  =  (1,  2,  3,  4). 

For  a  number  of  different  reasons  the  user  of  Subroutine  FFT  may 
want  to  alter  part  of  the  Fourier  frequency  spectrum  before  obtaining  the 
time  function  Xk.  Before  one  alters  the  frequency  spectrum,  however,  it  is 
important  to  know  that  the  second  half  of  the  frequency  spectrum  is  a  mirror 
image  of  the  first  half.  In  other  words,  the  value  of  Y^  in  the  first  half 
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equals  the  value  of  Yq_k  in  the  second  half,  where  q  equals  2n.  Therefore, 
if  some  of  the  values  of  are  changed  in  the  first  half  of  the  frequency 
spectrum,  the  corresponding  values  of  Y^  in  the  symmetric  half  must  similarly 
be  changed.  This  procedure  is  not  difficult,  but  it  is  necessary. 

REFERENCE 

Al.  E.A.  Robinson  Multichannel  Time  Series  Analysis  with  Digital 

Computer  Programs.  Holden-Day,  1967. 


C- . 

SUBROUTINE  FFT  (N,  X,  SIGN) 

C . - . . 

C  NMAX  =  LARGEST  VALUE  OF  N 

C  NONDUMMY  DIMENSION  M(NMAX) 

C  FOR  EXAMPLE,  IF  NMAX  =  8  THEN 

DIMENSION  M(8) 

C  DIMENSION  X(2**N) 

DIMENSION  X ( 2 ) 

COMPLEX  X,  WK,  HOLD,  Q 
LX  =  2**N 
DO  1  1=1, N 
1  M( I)  =  2**(N-I ) 

DO  4  L=1  ,N 
NBLOCK  =  2**(L-1 ) 

LBLOCK  =  LX/NBLOCK 
LBHALF  =  LBLOCK/2 
K  =  0 

DO  4  IBL0CK=1 ,  NBLOCK 
FK  =  K 
FLX  =  LX 

V  =  SIGN*6.2831853*FK/FLZ 
WK  =  CMPLX(COS(V) ,SIG(V)) 

I  START  =  LBL0CK*( IBLOCK  -  1) 
DO  2  1=1,  LBHALF 
J  =  ISTRAT  +  I 
JH  =  J  +  LBHALF 
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Q  =  X(JH)*WK 
X ( JH)  =  X(J)  -  Q 

2  CONTINUE 
DO  3  1-2, N 
II  =  I 

IF  (K.LT.M(I))  GO  TO  4 

3  K  =  K  -  M(I) 

4  K  =  K  +  M(  1 1 ) 

K  =  0 

DO  7  J  =  1  ,LX 
IF  (K.LT.J)  GO  TO  5 
HOLD  =  X(J) 

X(J)  =  X(K+1 ) 

X ( K+L )  =  HOLD 

5  DO  6  1=1 .N 
II  =  I 

IF  (K.LT.M(I))  GO  TO  7 

6  K  =  K  -  M(I) 

7  K  8  K  +  M(II) 

IF  (SIGN. LT. 0.0)  RETURN 
DO  8  1-1, LX 

8  X(I)  =  X( I)/FLX 
RETURN 

END 
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